External bias in the model of isolation of communities 
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We extend a model of community isolation in the d-dimensional lattice onto the case with an im- 
posed imbalance between birth rates of competing communities. We give analytical and numerical 
evidences that in the asymmetric two-specie model there exists a well defined value of the asym- 
metry parameter when the emergence of the isolated (blocked) subgroups is the fastest, i.e. the 
characteristic time t c is minimal. This critical value of the parameter depends only on the lattice 
dimensionality and is independent from the system size. Similar phenomenon was observed in the 
multi-specie case with a geometric distribution of the birth rates. We also show that blocked sub- 
groups in the multi-specie case are absent or very rare when either there is a strictly dominant specie 
that outnumbers the others or when there is a large diversity of species. The number of blocked 
species of different kinds decreases with the dimension of the multi-specie system. 

PACS numbers: 05.50.+q, 89.75.Hc, 02.50.-r 
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I. INTRODUCTION 

The question of imbalance or asymmetry is long known 
in many popular and fundamental phenomena in non- 
equilibrium statistical physics such as gradient-induced 
transfer of particles, heat or current The two-specie 
fermion mixtures and the impact of its population im- 
balance on the phase diagram of uniform superfluidity is 
also one of the multiple faces of the mentioned problem 
0,1- 

This issue emerges in a quite natural way in other 
fields: for example in ecology, where the imbalance can 
be related to inclusion of the non-native species in certain 
area [H, as well as to the absolute number of species in 
the branches of phylogenetic tress [5| or in general it can 
govern the whole evolution in a specific ecosystem Q . In 
economics the heterogeneity can be manifested for exam- 
ple by Zipf 's law in distribution of firm sizes 0, Q or in 
the statistics of order books [9(. 

Apart from the above mentioned quantitative sciences, 
there is also sociology, during the last two decades under 
frequent consideration of physicists which eventually led 
to the emergence of novel scientific entity — sociophysics 
10]. Its powerful tools have been widely used not only 
to face the problems of culture spreading [ll|, collective 
behavior of audienc e 113 , correspondence activity 13| or 
travel customs [H], Il5l| but very recently it also touches 
the issues of moral standards [16| and emotions 

EMU- 

By the same token, a significant attention is devoted 
to social balance. Its evolution on networks with both 
friendly and unfriendly relations has proven to provide 
signs of phase transitions [20( or complex energy land- 
scapes [21] . The complexity is inseparably binded with 
human nature with respect to the social imbalance: it 
is true that people are inequality- aware which is backed 
even by neural evidence [22] but then again they also tend 
to accept inequalities reflecting differences in individual 



achievements [231 ] . 

In this paper we draw attention to the issue of im- 
balance in a model of isolation of species (communities) 
extending the previously obtained results [2J] . Social iso- 
lation is a crucial problem which can be caused by various 
factors such as illness, imprisonment or emigration and is 
regarded to have a substantial impact on higher suicide 
rates [25j . In this context, isolation may also emerge as 
an effect of the social imbalance, just to mention large 
discrepancies between number of males and females in 
France after World War I (25|. Bearing that in mind, 
we shall try to show the impact that different methods 
of imbalance introduction may have on the number of 
isolated species. 

The paper's structure is organized as follows: in Sec. 
[IT] we describe in short the basic model of isolation, Sec. 
IIIII gives details of the ways the imbalance is introduced 
and Sec. II VI presents exhaustive analysis of the numerical 
simulations and the theoretical approach. The outcome is 
summarized in Sec. [V] along with some general remarks. 



II. BASIC MODEL 

Recently we introduced a simple model of community 
isolation (24j whose basic rules can be described in the 
following way: in each time step one puts a representa- 
tive of specie (f) or Q.) in a random, unoccupied node 
in a chain of N nodes. The specie can be an individual 
belonging to a given community or a person character- 
ized by a certain opinion. The probabilities of choosing 
either of the possible specie (community birth rates) are 
equal to 1/2. If n nodes filled with identical species (e.g., 
444-) § e f surrounded by individuals belonging to other 
community (e.g., t444t)i the nodes inside the cluster are 
called blocked and they no longer interact with the rest 
of the system. This model can be easily extended onto 



a case with multiple number of different species to > 2 
- then, similarly as in the two-specie case, the type of 
the specie is drawn from the uniform distribution (1, to). 
To form an isolated cluster, a set of identical species has 
to be surrounded by other identical species (e.g., 23332). 
The analytical treatment shows that the approximated 
equation for the number of blocked nodes Z for specific 
time of the simulation t is 



(to - l)t 3 
m 2 N 2 



(1) 



It follows, that regardless of the number of different 
species introduced in the chain, the number of blocked 
nodes grows roughly as t 3 . The form of Eq. (fTJ) gives the 
opportunity to directly obtain the characteristic time t c 
i.e., the time when the first isolated node appears 
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1/3 



N 2/3 . 



(2) 



The approach is quite consistent with the numerical sim- 
ulations and works also in case of other topologies such 
as lattices or random graphs [24j |. 



III. GENERAL ANALYTICAL DESCRIPTION 
OF THE ASYMMETRIC CASE 

In this paper we are investigating the asymmetric case 
of the model [2ij . Symmetry breaking is introduced as 
an external bias and is received by changing birth rates 
of different species occurring in the system, either by a 
simple imbalance in the two-specie case or by setting a 
specific specie probability distribution in the multi-specie 
one. We shall consider the model of isolation of commu- 
nities on a d-dimensional lattice. In this kind of graphs 
clusters of nodes sharing the same specie will be isolated 
when neighbors (in the sense of von Neumann neighbor- 
hood |26[) of all their elements are connected only with 
each other or with a different specie. An example of iso- 
lated cluster on the square lattice is presented in Fig. 

m 



A. Two-specie case 

First, we consider two species (f and ],) which we put 
into d-dimensional hypercubic lattice with a total number 
of sites N = L d . Occurrence probabilities p-f and p± of 
both species are non-equal and given by 



1 

Pt= 2 



1 

PI =2 



(3) 



where e G (0,1/2) is the symmetry breaking param- 
eter. These probabilities are microscopic parameters, 
which describe model evolution. Starting from a lattice 
of empty nodes, after t time steps we have the following 
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FIG. 1: (Color online) An example of cluster of size 3 on 
square lattice (d — 2) with N = 20 nodes and two different 
species (light and dark filled circles), after t — 16 time steps. 
Nodes marked with crosses are blocked. 



probabilities that a randomly picked node is occupied 
with the specie t or 4- 

Prob(t) = ^Pt, Prob(|) = ^. (4) 

Now let us consider numbers of blocked nodes , of 
both species at time t. These numbers are well approx- 
imated by numbers of individual blocked nodes, i.e. by 
numbers of blocked clusters of size one of both species. 
When the total density of all blocked nodes is small 
{p + Z± < t) then 



Z^ 



Z\ = (L- 2) d Prob 2d (|)Prob(t), (5) 



Z l Z[ = (L - 2) d Prob 2d (t)Prob(|). 



(6) 



Substituting the values Prob(f) and Prob(4,) with Eq. 
(j4|) and taking into account ((3]) the above relations for 
L 3> 2 may be expressed as 



t 



2(2+1 



z± 



N 2d 



i2d+l 



N 2d 



1 



2d-l 



2d-l 



(7) 



(8) 



The above equations allow us to calculate characteristic 
times t\ and tj; when the first representative of either of 
both species emerges. Putting Z* = 1 (analogously for 
4-) we obtain after a short algebra 



N 



2,1 



N 2(2-1 



N 



2d 



(\-c 2 )(h + c) 



\ 2(2-1 



(9) 



(10) 
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FIG. 2: (Color online) (a-b) Number of blocked nodes Z-f and Z^ versus symmetry breaking parameter e on a chain (a) and 
on a square lattice (b); circles (Zf) and triangles (Zj.) correspond to simulation results and lines are calculated using Eqs 
(J7]|8j) . (a) Simulations performed for t = 7V/4 with N = 50000 nodes, averaged over 100000 runs, (b) Simulations performed 
for t — N/2 with N = 10000 nodes, averaged over 1000 runs. The regions on left-hand side of the vertical line of each plot 
correspond to positive values of the derivative dZ^/de while in those on the right-hand side dZ^/de < 0. The line dividing two 
regions is drawn for e = e* = 1/6 for the one-dimensional case and e = e* =3/10 for the two-dimensional as predicted by Eq. 



Q5) . (c-d) Log-linear plots of the characteristic times tc (circles) and t 1 ^' (squares) for d = 1 (c) and d — 2 (d). Data points 
are taken from numerical simulations while solid lines come from Eqs (|51 ll0[l , (c) N = 50000, averaged over 100000 runs, (d) 
N = 10000, averaged over 10000 runs. The insets in both plots show in detail the range for which tc takes the minimum value 
(e = e* = 1/6 for one-dimensional case and e = e* = 3/10 for the two-dimensional) marked by vertical dotted lines. 
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B. Multi-specie case 



tain that 



If we consider more than two species, we can use prob- 
ability distribution Prob(i) instead of the previously in- 
troduced probabilities Prob(f) and Prob(4_). It requires 
changing microscopic probabilities p^ and p± with the 
value of pi. The connection between these two sets of 
variables is given by the equation 



Prob(z) = — 



(11) 



Z w w Zf w N J2 Prob(fc) 2d Prob( i ) (12) 

k=i, k=£i 

for i = 1,2,..., which, taking into account (fTTj) becomes 



Pi\Epl d -Pi d 



(13) 



Putting Z« = 1 one obtains a set of characteristic times 
for each specie i 



for i = 1, 2, ... (the case pi independent of i corresponds 
to the symmetric problem considered in (24|). Using the 
same approximation as in the two-specie case we can ob- 



N 



2<l 



Pi^pf-Pf) 



i = l,2,. 



(14) 
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FIG. 3: (Color online) Log-log plot of the global characteris- 
tic times t c versus parameter q of the geometric distribution 
USJ for Ni = 15625 (empty symbols) and N 2 = 65536 (rilled 
symbols). Data points are taken from numerical simulations: 
squares are d = 1, circles are d — 2, triangles are d — 3 and 
diamonds are d — 4; all data points are averaged over 1000 
runs. Solid lines come from Eq. (|17[1 and the horizontal solid 
line is drawn for t c = Ni whereas the horizontal dotted one for 
t c = Nn. The vertical solid line marks q m in ~ 1/Ni while the 
vertical dotted line is drawn for qw„ = 1 /N% (see description 
in text). 



IV. NUMERICAL SIMULATIONS AND 
DISCUSSION 

Here we are presenting the comparison between the an- 
alytical approach given in the previous section and the re- 
sults obtained from the numerical simulations performed 
for topologies of <i-dimensional lattices. 



A. Two-specie case 

While comparing Eq. (JTJ for m — 2 and Eqs. (05]) for 
specific value of e one can notice that the introduction of 
the symmetry breaking does not have any effect on the 
dependence of the number of blocked nodes with regard 
to the time of the evolution. However, investigation of 
and as a function of the parameter e for a fixed value 
of time t brings some not so obvious results. Figure [5^-b 
shows the dependence of the number of blocked nodes 
of each specie versus the symmetry breaking parameter 
e obtained in the numerical simulations compared with 
the theoretical expectations given by Eqs (05J) . The dis- 
crepancies between simulation results and the theoretical 
approach can be justified by the approximations used in 
obtaining Z l and Z f . Still, one can immediately spot 
the main difference between those two quantities: Z' is 
monotonic while Z^ first increases, reaches a prominent 
and well defined maximum and then drops down. This 



observation is backed with a simple analysis of Eqs. (J7J- 
|5J| : in fact the derivative dZf /de < for the whole range 
e G (0,1/2) and in the case of Z^ there is a maximum 
value for 



1 2d - 1 

2 2d + 1 ' 



(15) 



The decrease of Z^ is rather obvious as it is the dominant 
specie according to Eq. (|3J) - the higher is the number 
of its individuals introduced in the system the smaller 
is the probability of this specie being blocked. On the 
other hand as it concerns Z^, Eq. (|15p suggests that 
there exists a specific value of e for which the number of 
blocked individuals reaches the highest rate. It follows 
that this value depends only on the dimensionality of the 
system (i.e., the number of the nearest neighbors). This 
phenomenon is even more pronounced while examining 
the characteristic times tj and i* versus e presented in 
Fig. [5J>d. The first quantity exhibits a constant growth 
and the second one possesses a clear minimum for e = 
e* which is consistent with the maximum values for Z^ 
observed in Fig. [2^,-b. 



B. Multi-specie case 

For the asymmetric multi-specie case we used the geo- 
metric probability distribution 



Pi = 1 



1,2,. 



(16) 



This very case of probability distribution has been cho- 
sen as an example as it quite easily yields the analyti- 
cal approach. Moreover it can be transformed directly 
into continuous exponential distribution assuming that 
q = exp(— a). First we shall discuss the issue of the 
global characteristic time t c . It can be approximately 
calculated by assuming that the total number of blocked 
species at time t is equal to Z(t) = = 1. Using 

Eq. (fTB"]) the global characteristic time can be expressed 
as 



tr. = 



N 



2<l 



(l_g)2<J+l 



(17) 



Figure [3] presents the characteristic time versus parame- 
ter q given for different values of d and N. In each case 
the curve has a similar shape, exceeding t c = N for both 
small and large values of q (horizontal lines in Fig. UJ) 
with a well defined minimum between. Moreover, for a 
specific value of N the curves, regardless of the dimen- 
sionality of the system, seem to intersect in one point. 
In fact, a closer analysis reveals that for g<l Eq. (fT7|) 
takes a form of t c = (N 2d / , g) 1 /( 2rf + 1 ) resulting in the in- 
tersection point q m in = l/N, shown as vertical lines in 
Fig. UJ A heuristic explanation of this fact can be also 
derived in the following way: as long as the probability of 
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FIG. 4: (Color online) (a-b) Log-log plots of characteristic times tc versus parameter q of the geometric distribution (|16[) for 
d = 1 (a) and d — 2 (b). Data points are taken from numerical simulations: squares represent specie i = 1, circles - i = 2, 
upward triangles - i = 3, downward triangles i = 4, diamonds - i = 4 and stars - i = 6; solid lines come from Eq. (|19l) . 
Theoretical curves for i = 1 and i = 3 in one-dimensional case as well as i = 1 and i — 5 in two-dimensional overlap. In both 
cases N = 15625 and the horizontal dotted line is drawn for ti^ — N. Numerical data are averaged over 1000 runs, (c) Log-log 
plot of the characteristic time ti^ calculated from Eq. (|20[) for q = q* against specie number i. Squares are d = 1, circles are 
d = 2, triangles are d — 3 and diamonds are d = 4; in each case JV = 15625. Solid lines are guidance to eye and the horizontal 
dotted line is drawn for t c = N. 



drawing the second specie (i — 2) is above 1/N there is a 
statistical chance of it appearing in the system and thus 
being blocked by the overwhelming first specie (i = 1). 
As soon as the probability drops below that value there 
is only one specie and so the blocking is impossible. On 
the other hand, when q approaches 1, the number of dif- 
ferent species is relatively high and all the probabilities 
Pi are close to p\ . After crossing a certain value q max the 
number of the individuals of the first specie is too low 
to make possible any blocking on the lattice. The obvi- 
ous necessary condition preventing this scenario is that 
Np 2d (p 2 +P3 +Pa + ...) = 1 which leads to Np 2d m 1 and 
eventually gives q m ax ~ 1 — A r ~ 1 / 2d . Finally Fig. |3] sug- 
gests that there is some specific value of q for which the 
characteristic time is the lowest; assuming that q 2d < 1 
one can estimate this value with 



1 



q = 



2d+l 



(18) 



We stress here that the critical value of q* given by Eq. 
(fT5|) is independent on the system size N similarly as the 
critical value of the e* parameter given by Eq. (1151) . 

One can also focus directly on the question of char- 
acteristic times ic of various species. Substituting Eq. 
fll~4l with (|16l) one gets the general formula for the charac- 
teristic time of i-th specie being blocked in the topology 
of the d-dimensional lattice 



N 2d (l - q Id ) 



L (l - g )2d+l (1 - g2d(i-l) + q 2di^ 



(19) 

(i) 

The comparison of the characteristic times t c versus pa- 
rameter q obtained in numerical simulations and the the- 
oretical expectations given by Eq. (|19l) is presented in 



Fig. |4J The plots for tc bear close resemblance to those 
obtained for t c . Following a similar line of thought as in 
the case of the global characteristic time it is possible to 
estimate the crucial points of these curves. First of all, 
Eq. (fT9"| can be approximated as 



(1) ^ (N 2d q- 2d ) 2d+l 
c ~ T^q 

(l ) _ (Af 2 V-')WT 
c ~ 1-q 



(20) 



1 = 2,3,... 



which is valid for all q G (0,1). However, in order 
to estimate q <^ 1 when the above function intersects 
with # = N, one uses a further approximation i.e., 
t 0O = N 2d q (i-i)/(2d+i)_ It c ff ccts i n qmte straightforward 
formula for the intersection point being q m i n — N 1 '*- 1- *', 
which suggests that the intersection point should be de- 
pendent only on the number of nodes in the lattice, while 
its dimension d is irrelevant. Figure 4 also gives clear ev- 
idence that the characteristic times for the first specie 
behave differently than those for i > 2, in both cases it is 
the specie i = 2 which gets blocked first for smaller values 
of q. Moreover in the case of d = 1 the results from t^ 
cover with i^ 3 ' and in the case of d — 2 with t c 5 ^ ■ This last 
fact is fully comprehensible after comparing both equa- 
tions in (|2TJ)) where, after short algebra, one gets that the 
characteristic time for the first and (2d+ l)-th specie are 
the same. One can also notice comparing Eqs (|17j) and 
(fl"9f that it is in fact t c that plays the dominant role 
in creating the shape of the global characteristic time t c 
and may be used as its good approximation. 

Taking into account Eq. (|20)) it is possible to esti- 
mate the maximal specie number that gets blocked in 
the environment for a fixed value of the system size N . 
The curves presented in Fig. 0Ji-b indicate that for each i 
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there is a specific value q* for which the function takes 
its minimum. It follows, that as long as tc (q*) < N the 
specie will be blocked, at least for q — q*. Closer analysis 
of Eq. (J20]) leads to 



n * _ 2d 
1l - 4d+l 

a* — 1=1 

It i+2d 



i = 2,3, 



(21) 



After substituting Eq. (|20l) with the above values of q*, 
one gets the value of characteristic time in the minimum. 
A corresponding plot for different values of dimension d 
is shown in Fig. HJ;. It gives immediately the idea of 
the fast restriction of the blocked specie number with 
system's dimension: while for d = 1 the value of i can 
be as large as 70, in the case of d = 4 it substantially 
drops down to 6. Furthermore, Fig. 0t underlines again 
the specific role of the first specie showed in the previous 
paragraph. 



trivial results. We have found that both in the two-specie 
case where a parameter e governs the symmetry break- 
ing and in the multi-specie case where the numbers of 
each entity are given by the geometric distribution there 
exists some specific and well defined value of the con- 
trol parameter giving a minimum of the characteristic 
time t c . In general the value of this parameter is depen- 
dent only on the dimensionality of the lattice. We have 
also shown that the requirement of a non-blocked system 
in the multi-specie case leads to two, somehow opposite 
conditions: either there has to be a strictly dominant 
specie, outnumbering the others or the diversity should 
be very large. In the end we have also given evidence 
that the number of blocked species of different kinds de- 
creases with the dimension of the system. The presented 
results can be easily further generalized to the cases of 
other topologies (e.g. complex networks) and other kinds 
of biases. 



V. CONCLUSIONS 

In this paper we have extended the simple model of 
community isolation onto the case with the symmetry 
breaking. Our calculation and numerical simulations 
show that even a simple way of introducing the exter- 
nal bias between species can lead to interesting and non- 
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